Cohesion and conductance of disordered metallic point contacts 



J. Burki 

Institut de Physique Theorique, Universite de Fribourg, CH-1700 Fribourg, Switzerland and 
Institut Romand de Recherche Numerique en Physique des Materiaux, CH-1015 Lausanne, Switzerland 

C. A. Stafford 

Physics Department, University of Arizona, 1118 E. 4th Street, Tucson, AZ 85721 



X. Zotos 

Institut Romand de Recherche Numerique en Physique di 



Materiaux, CH-1015 Lausanne, Switzerland 



D. Baeriswyl 

Institut de Physique Theorique, Universite de Fribourg, CH-1700 Fribourg, Switzerland 
(Submitted to Phys. Rev. B, 25 February 1999, resubmitted 4 May 1999) 

The cohesion and conductance of a point contact in a two-dimensional metallic nanowire are in- 
vestigated in an independent-electron model with hard-wall boundary conditions. All properties 
of the nanowire are related to the Green's function of the electronic scattering problem, which is 
solved exactly via a modified recursive Green's function algorithm. Our results confirm the valid- 
ity of a previous approach based on the WKB approximation for a long constriction, but find an 
enhancement of cohesion for shorter constrictions. Surprisingly, the cohesion persists even after the 
last conductance channel has been closed. For disordered nanowires, a statistical analysis yields 
well-defined peaks in the conductance histograms even when individual conductance traces do not 
show well-defined plateaus. The shifts of the peaks below integer multiples of 2e 2 //i, as well as the 
peak heights and widths, are found to be in excellent agreement with predictions based on random 
matrix theory, and are similar to those observed experimentally. Thus abrupt changes in the wire 
geometry are not necessary for reproducing the observed conductance histograms. The effect of 
disorder on cohesion is found to be quite strong and very sensitive to the particular configuration 
of impurities at the center of the constriction. 

PACS numbers: 73.23.-b, 73.20.Dx, 72.80.Ng, 73.40. Jn 



I. INTRODUCTION 

In a pioneering experiment published in 1996, Rubio, 
Agra'it and VieiraEl simultaneously measured the force 
and conductance during deformation and rupture of gold 
nanocontacts. They observed nano-Newton force oscilla- 
tions correlated with conductance steps of order 2e 2 /h. 
Similar results were obtained independently by Stalder 
and Diirig.u An explanation of these force fluctuations, 
based on the response of the conduction electrons to the 
mechanical deformation of the contact, was proposed by 
Stafford et al.B In this paper, we give an exact numerical 
solution to the model of Ref. ||, allowing the treatment 
of any shape of the constriction and the inclusion of dis- 
order. 

Conductance quantization has been widely.— stud- 
ied in the nasi .decade, both experimentallyQtll and 
theoretically.Ellljt3 On the experimental side, it was first 
measured in a two-dimensional electron gas split-gate and 
was understood to be a consequence of the quantization 
of the transverse motion.oO Each transverse energy sub- 
band defines a conduction channel that contributes 2e 2 /h 
to the conductance. More recently, conductance quan- 
tization in units of 2e 2 /h has been observed in three- 



dimensional metallic cmantum wires using various exper- 
imental techniquesJjEj In this case, the quantization is 
not so clear-cut, and measurements are in general less re- 
producible. Therefore, a statistical analysis of a large set 
of experimental data is usually made. The resulting con- 
ductance histograms show well-defined peaks with cen- 
ters shifted below integer multiples of 2e 2 /h. The shifts 
can be corrected by subtracting a series resistance of a 
few hundred Ohms. This resistance was originally at- 
tributed to the bulk contacts, which are noUuaxt. of the 
nanocontact and should thus be subtracted. I'LTtiJ-it has 
recently been argued, based on theoretical work,Eil that 
this resistance could also be caused by disorder in the 
nanowire, in agreeipaHt. with previous studies of quan- 
tum point contacts. EJJEJ 

On the theoretical side, the problem of conductance 
quantization has been investigated using anyariety of 
techniques, both in two and three dimensions BEd^a It is 
found that conductance plateaus are well defined only in 
constrictions^tliat| are smooth on the scale of the Fermi 
wavelength.li3ii3'tj Furthermore, the plateaus can also 
be affected by thermal fluctuations or disorder. Thex=. 
mal effects, which smoothen out the conductance steps,L3 
are important in the case of a two-dimensional elec- 
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tron gas ,00 where the spacing between subbands is of 
the order of lJsTmrbwti negligible in the case of metallic 
point contactSjIila-Ertil where the spacing between sub- 
bands is of the order of 10 A K. The effect of disor- 
der on condiijct.a.nce quantization hasJacen studied both 
analyticallyBBB and numerically.lOEj Disorder was 
found to reduce the conductance compared to the ballis- 
tic case, leading to shifts of the peaks in the conductance 
histograms similar to those observed experimentally. We 
note that most of the theoretical conductance curves were 
obtained by changing the Fermi energy, in contrast to the 
experimental situation for a metallic nanowire, where the 
Fermi energy is an invariant property of the material and 
the shape of the contact is modified. 

The structural transformations of metallic nanocoB-. 
tacts were first studied by Landman and coworkers,cJ 
and later by othersJ3llj using molecular dynamics. These 
simulations suggest that the elongation of a connec- 
tive neck proceeds through a sequence of abrupt struc- 
tural transformations, involving a succession of elastic 
and yielding stages. Within this approach, the con- 
ductance is expected to change abruptly due to atomic 
rearrangements,Bll3 and. the force is expected to follow 
a sawtooth behaviorJao The molecular dynamics-simu- 
lations give rise to the following interpretationErc3l23 of 
experiments: During elastic stages, the conductance is 
essentially constant at a quantized value and the magni- 
tude of the force increases with a constant slope, while 
the yielding stages consist of abrupt relaxations of the 
atomic structure, giving rise to a sudden change of con- 
ductance and force. It was asserted that both the elastic 
and transport properties of narrow metallic constrictions 
can be understood on the basis of atomic rearrangements, 
calculated in the framework of classical lattice dynamics. 

It thus came as a surprise when it was discovereda that 
a simple free electron model, which neglects the discrete 
atomic structure of the constriction, is able to reproduce 
both the oscillations in the tensile force and the conduc- 
tance plateaus of Ref. [|. In this model, each conduc- 
tance channel is viewed as a long chemical bond which 
is stretched and broken, giving rise to the observed force 
oscillations—Similar results were obtained by van Ruitcn- 
beek et alE3 using a f-eee-electron model and by Yan- 
nouleas and Landman,t3 using a slightly more realistic 
jellium model. However, these studies dealt with closed 
systems, while the experimental nanowire is a contact 
between two macroscopic pieces of metal, and is thus an 
open system. Since mesoscopic effects can be very differ- 
ent in the canonical and grand canonical ensembles, this 
is an important distinction. 

In this paper, we give an exact numerical solution of 
the free-electron model proposed in Ref. ^. Our formal- 
ism allows us to extend the analytical results of Ref. || 
by treating both non-adiabatic constrictions and disor- 
dered wires. We find that the cohesion of short (non- 
adiabatic) constrictions is increased compared to that of 
long constrictions. For disordered systems, we obtain 
conductance histograms which look very similar to those 



observed experimentally. Our results also suggest that 
the cohesion of nanowires is quite sensitive to disorder. 

The rest of this paper is organized as follows: In sec- 
tion [n| we describe our model in which conductance and 
cohesion are interconnected through the scattering ma- 
trix. Section |l| presents the numerical method, while 
the results are discussed in section IV. A short summary 
is given in section 



II. MODEL AND BASIC QUANTITIES 

This section introduces our free-electron model. Its 
simplicity allows to treat both the conductance and the 
metallic cohesion in a unified way. 



A. Assumptions 

Our model is based on several simplifying assumptions. 
(i) Free electrons: Gold is known to be rather well de- 
scribed by the model of pee electrons where the atomic 
potentials are neglected.E^I This indicates that the de- 
tailed configuration of atoms has little effect on both elec- 
tronic and cohesive properties. By adopting the free elec- 
tron model, we eliminate the kinetics of atomic rearrange- 
ments, which seem to play a role, in actual experiments, 
at least for thick constrictions.E3 (ii) Independent elec- 
trons: We neglect electron-electron interactions. Recent 
calculations, one based on an extended Thomas Jiermi 
scheme,tll the other on the Hartree approximation,^ sug- 
gest that interactions have little effect on the quanti- 
ties considered here. However, more careful studies will 
be needed to substantiate this conclusion, as interac- 
tion effects arc particularly subtle in one-dimensional 
systems and thin wires, (hi) Smooth constriction^ Ac- 
cording to molecular dynamics simulationsulJt§E3 and 
experimentsj23'Ej the shape of elongated gold nanowires 
is quite regular. Therefore we model the contact as a 
smooth constriction, which changes continuously upon 
deformation. (iv) Constant volume: The volume of 
the constricted part is assumed to be conserved during 
the deformation- in agreement with_molecular dynam- 
ics simulations ril Recent experimentsE3 indicate that this 
assumption is not always valid, for instance in the case 
of mono-atomic chains, (v) Hard-wall boundary condi- 
tions: Electrons are confined to a wire due to a smooth 
attractive ionic potential. We approximate the potential 
in terms of hard-wall boundary conditions, keeping in 
mind that the radius of the boundary will be larger than, 
the effective radius defined by the electronic density.EZl 
(vi) Zero temperature: The thermal population of elec- 
tronic subbands above the Fermi energy is negligible for 
nanowires. Therefore we consider the zero-temperature 
limit, where inelastic scattering processes, e.g. due to 
phonons or electron-electron collisions, are absent. 
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B. Model 

The physics of the problem is similar for two- 
and three-dimensional wires except that, for a two- 
dimensional system, (i) the transverse energy levels are 
nondegenerate, and therefore the steps in the conduc- 
tance curves are all of height 2e 2 /h, and (ii) the sur- 
face energy is reduced as compared to that of a three- 
dimensional wire; this changes the overall slope of the 
force curves. We restrict ourselves to two dimensions, 
where the computational effort is much lower. 

Our model wire is described by the eigenvalue equation 



2m 



(d 2 x + d 2 )+V dis (x,z) 



ip(x, z) = Eip(x, z), (1) 



where x and z are the transverse and longitudinal co- 
ordinates, respectively, and Vdi 8 {x,z) is a potential due 
to disorder. We use the boundary condition -0 = for 
x 2 = r 2 (z), where r(z) defines the geometry (see Fig. 
pi). In the initial configuration, there is no constriction 
Jr(z) = R) and the disorder is restricted to a finite part 
of the wire (namely to \z\ < (Lq/2 + L^is), where Lq is 
the initial length of the deformable part). After elonga- 
tion by AL = L — Lq , the shape of the deformable part 
is chosen to be, for \z\ < L/2, 



r(z) = Rmin + {R~ Rmin) (3u 2 - 3u 4 + U 6 ) 



(2) 



where u — | =jr- 1 and 2R m i n is the width of the narrowest 
part of the constriction, related to L by the condition 
of constant volume. The specific choice of the function 
r(z) is not crucial,u as has been checked by using other 
shapes. 





-dis 



L=L„+ AL 



are moved in such a way that their local concentration 
remains constant. Since we are dealing with an open sys- 
tem, we use the grand-canonical ensemble and take the 
chemical potential to be the Fermi energy of the macro- 
scopic wires connected to the nanocontact. Throughout 
this paper, we consider the case of gold, with a Fermi 
energy ep = 5.5 eV. 



C. S matrix 

The formalism developed by Stafford et al.i allows one 
to treat both transport and mechanical properties of open 
mesoscopic systems on an equal footing, in terms of the 
electronic scattering matrix S. This formalism is general 
and can be applied to any non-interacting open system 
connected to an arbitrary number of leads. 

The conductance for the present two lead system is 
calculated using the Landauer formulaEj 
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(3) 



where t is the transmission matrix at the Fermi energy 
(directly related to the S matrix). 

The force is obtained from the variation of the appro- 
priate free energy, the grand canonical potential Q(ep, L), 
with respect to the length L of the deformed region. For 
a non- interacting system of electrons, one has 



n(e F ,L)=J dE{E -e F )D(E), 



(4) 



where the density of states D(E) is related to the S ma- 
trix byBH 



(5) 



-dis 



The S matrix provides a fundamental link between con- 
ductance and cohesion, and it is amenable to analytical 
calculations B S is directly-related to the retarded elec- 
tronic Green's function,E3t3 which may be calculated nu- 
merically, as explained in the next section. A closely re- 
lated formulation of the free energy in. terms of the S 
matrix was given by Akkermans et al.E3 



FIG. 1. Schematic diagram of a constriction in a 2D 
quantum wire. Electrons are confined along the z axis by 
a hard- wall potential at x = ±r(z). Localized impurities of 
potential V(x — Xj,z — Zj) are distributed randomly in the 
constriction and within a part of length Ldis on each side of 
it. 



The potential V dls (x,z) = £\ V(a 



h 3 ' ' 



;) de- 



scribes randomly distributed impurities, with a given 
density ri$. It can represent structural defects or real 
impurities. Similar effects are expected to arise due to 
surface roughnessH During deformation, the impurities 



III. EXTENDED RECURSIVE GREEN'S 
FUNCTION METHOD 

In this section, we describe our numerical procedure. 
It consists of computing recursively the Green's function, 
from which both the transmission coefficients and the 
density of states are obtained. 
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A. Green's function for a wire of constant width 

In order to illustrate our method, we first consider 
an infinite, two-dimensional wire of constant width 2R. 
In order to proceed numerically, the original continuum 
model is replaced by a discrete lattice, with a width of M 
sites, while the disordered part is restricted to a length 
of N sites. This finite region consists of slices, numbered 
from 1 to N, each slice representing a cross section with 
M sites. The discretized Hamiltonian (for the disordered 
part of the wire) is then a (NM x NM) matrix, with 
M x M sub-matrices for the i-tb slice and M x M 
sub-matrices for the hopping terms between neigh- 
boring slices i and j. 

The Recursive Greeajs Function (RGF) method, devel- 
oped by MacKinnon ,l2I is based on the Dyson equation 
for the Green's function. At each step of the iteration, 
only the relevant parts of the Green's function are re- 
tained. Thus the RGF method allows one to calculate 
both the density of states and the transmission coeffi- 
cients at low memory cost. Applied to our problem, the 
method allows us to construct the Green's function of 
the disordered part of the wire slice by slice. The infinite 
ordered parts are included through boundary conditions 
for the first and last slices of the disordered region (see 
Ref. [h] for details) . The advantage of this method is that 
one only needs to keep track of a few M x M matrices 
instead of the (NM x NM) matrix representing the full 
Green's function. 



B. Generalization to a wire of variable width 

In order to generalize the RGF method to the case of a 
wire with varying thickness, we first apply a scale trans- 
formation. Consider the eigenvalue equation (Qj) together 
with the boundary condition tp(x,z) = for |ac| = r(z). 
The change of coordinates 



R 

'r(zY 



(6) 



transforms the wire geometry into a stripe of constant 
width 2R, at the cost of a more complicated differential 
operator, 



8l = \ 



R 2 



2 {r 2 (z) 



x 2 p 2 ,dl 



dj - {xp,d s d s } - 
3 



1 dp 

2 dl' 



(7) 



where p(z) = ■4=log[r(z)] and {.,.} is the anticommuta- 
tor. This operator is well-behaved if r(z) is twice contin- 
uously differentiable. 

The stripe geometry now allows a straightforward dis- 
cretization. For the case of interest, where r(z) = R 
outside of the constriction, the RGF method can again 
be used, as described above. The transformation (0) has 



the additional advantage that the narrowest part of the. 
wire, where the main effects considered here originate, a 
is scanned with the finest grid. 



C. Calculation of conductance and force 



Once the retarded Green's function Gy is determined, 
the transmission matrix in the Landauer-Buttiker for- 
mula [Eq. (||)] can be easily calculated (see Ref. |d|). The 
resulting equation for the conductance is found to be 



dEu. dE v 



dk dk 



(Gol) 



(8) 



where E v (k) is the dispersion relation for the v-th sub- 
band. The Green's function also yields directly the den- 
sity of states 



D(E) = ImTrG(E), 

7T 



(9) 



from which the grand canonical potential [Eq. (0)] is ob- 
tained by integration. The tensile force is then calculated 
by differentiating the energy f2 with respect to the length 
L of the constriction, 



F(L) = 



dCl(L) 
dL ' 



(10) 



For the conductance, it is sufficient to know the Green's 
function at the Fermi energy. The evaluation of the force 
requires much more work, as all occupied electronic states 
contribute to the grand canonical potential [Eq. (Q)] . It is 
thus necessary to recalculate the Green's function many 
times, and this is the part that uses most of the compu- 
tational effort. Fortunately, the energy Q.(L) turns out 
to be a smooth function of L, so that the numerical dif- 
ferentiation can be carried out for rather large steps AL. 



IV. RESULTS 

In Ref. ||, we solved the free-electron model with- 
out disorder analytically, using both the adiabatic and 
WKB approximations. In this section, the validity of 
our analytical results is confirmed numerically for long 
and clean constrictions. We then proceed to the case of 
short necks, where the adiabatic approximation breaks 
down. Subsequently, impurities are introduced in terms 
of randomly distributed short-ranged potentials. Both 
non-adiabaticity and disorder are found to produce in- 
teresting new effects. 



A. Comparison with analytical results 

We first want to check the analytical results of Ref. 
|3|. We thus consider a clean wire with a long constric- 
tion, i.e., with a small enough ratio R/Lq, where 2R and 
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Lq are, respectively, the width of the wire and the initial 
length of the constriction. The width is fixed to UfR = 6. 
We have varied the lattice constant a of our mesh in or- 
der to study possible discreteness effects. We have found 
that these effects are negligible for kpa < 0.7. All re- 
sults presented below were obtained for lattice constants 
satisfying these conditions. 

In Ref. ||, it was shown that within the adiabatic and 
WKB approximations the force is invariant with respect 
to a stretching of the geometry r{z) — > r(Xz), and thus 




L/L 



FIG. 2. Comparison of analytical (WKB) and numerical 
(RGF) results for a wire of width 21cfR = 12. For the conduc- 
tance, the WKB result (/cf£o = 120, solid line) is compared 
with two RGF results (JcfLq = 120, dotted line; IcfLq — 12, 
dashed line). For the force the WKB result (solid line, inde- 
pendent of Lq) is displayed together with several RGF results 
{k F L = 120, 36, 24, 12, from top to bottom). 



is independent of the initial length Lq. We therefore com- 
pare in Fig. H the analytical results for a single initial 
length, kpLo = 120, with numerical results for different 
Lo- For kpLo = 120 the conductance curves are almost 
indistinguishable, and the numerical result for the force 
comes very close to the analytically calculated curve. The 
analytical approach is thus seen to be justified for a long 
constriction. The behavior for shorter constrictions will 
be discussed in the next subsection. 

An intriguing effect is observed for an elongation 
L/Lq > 1.5. The force, expected to change sum as 
soon as the last conductance channel is closed,B i.e., 
for L/Lq ps 1.5, is found to remain attractive until 
L/L w 1.6. For the case of gold, this would correspond 
to a stretching of about 5A beyond the point where the 
conductance becomes vanishingly small n This surprising 
effect may actually have been observed.EJ 

We have checked whether this persistent force origi- 
nates from our arbitrary choice of the shape of the con- 
striction. Therefore, we have determined the shape that 
minimizes the macroscopic part of the free energjo for a 
fixed elongation. With our constraint, this is equivalent 
to minimizing the boundary length of the constriction, 
subject to the condition of constant area. We find that 
in this case the effect is even stronger than in other ge- 
ometries. Nevertheless, in three dimensional wires, this 
persistent force depends strongly on the shape and some- 
times even disappears. Further work is needed to clarify 
the origin of this interesting effect. It is worthwhile to 
mention that for this optimized shape the conductance 
steps turn out to be more pronounced and closer to ex- 
periments than for other shapes. 





1.6 ^ 



1.6 ^ 



L/L 



L/Ln 



FIG. 3. Conductance (a) in units of 2e 2 /h and force (b) in units of ef/^f as functions of elongation L/Lo and opening angle 
for a fixed width 2kpR = 12. 
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B. Non-adiabatic constrictions 

Our next goal is to go beyond the limitations of the adi- 
abatic approximation and to study short constrictions. 
To this end, we consider clean systems with a fixed radius 
kpR = 6, and vary the initial length Lq. We characterize 
the geometry of a wire by its opening angle 9, defined by 
tan0 = 2R/L . 

In the adiabatic approximation the wavefunction of Eq. 
(|l|) is factorized, ip(x, z) — <jf(x)x(z)> and terms contain- 
ing derivatives of 4> z {x) with respect to z are neglected. 
This solution corresponds to a set of decoupled chan- 
nels. In the present context, the approximation is valid 
if \-^r\ <C k, where k is the wave vector of the incident 

I r dz I 1 

electron. For the case of the conductance, where only 
states at the Fermi energy are involved, this condition is 
equivalent to 9 <C tt/2. Figure ||(a) shows the conduc- 
tance versus elongation L/ Lq and opening angle 9. These 
results for the. conductance are consistent with those of 
Torres et al.£3 for an adiabatic wire (9 <C tt/2), the con- 
ductance is well quantized in units of 2e 2 /h, with sharp 
steps between the plateaus, while for non-adiabatic wires 
the edges are rounded off and the plateaus are no longer 
well defined. The effects of the geometry on the force are 
shown in Fig. ||(b). One notices that as 9 increases, the 
average force becomes larger and the force oscillations 
are slightly suppressed. The overall increase of the force 
is due to the increased effect of surface tension in short 
constrictions. In the two-dimensional case, the surface is 
one-dimensional, and-jthe surface energy is proportional 
to the circumferenceo 

fisurf = ^fdA = [ L dz [1 + [dr/dzf] 1/2 . (11) 
6af oaf Jo 

In the WKB approximation, the circumference is approx- 
imated by dAwKB = 2L, which is a good approximation 
provided (dr/dz) 2 <C 1, as pointed out in Ref. |[ How- 
ever, this approximation leads to an underestimate of 
the surface tension in short constrictions, where \dr/dz\ 
is large. The increased effect of surface tension in short 
constrictions with a special wide-narrow-wicLe geometry 
was recently discussed by Kassubek et al.Ej As to the 
force oscillations, we note that they remain well-defined 
provided that one can resolve the conductance plateaus. 
Thus we attribute their decrease with increasing values 
of 9 to enhanced tunneling. 

We remark that wires with larger 9 can be stretched 
more than ones with small 9, that is to say, the cohesive 
force remains negative (attractive) further past the point 
where the conductance goes to zero. Thus the remanent 
cohesion is enhanced for non-adiabatic constrictions. 



C. Effects of disorder 

We now turn to the study of disorder effects. We first 
calculate the conductance for a single impurity configu- 



ration, then construct histograms from several hundred 
samples and compare them both with experimental re- 
sults and previous numerical and theoretical work. Fi- 
nally, we study the effect of disorder on the tensile force. 

1. Conductance for a single impurity configuration 

We model an impurity at rj = (xj,Zj) by a local po- 
tential V(x, z) = jS(x — Xj)8(z — Zj). The disorder for a 
density n.; of impurities can be characterized by the mean 
free path in the Born approximation, £ = ^2^-5 • Fig. |] 
shows conductance curves for a disordered wire with 7 oc- 
cupied channels (kpR = 12, kpLo — 60), with an impu- 
rity concentration Xprii = 0.27 both within the constric- 
tion and over a length kpLdis = 36 on each side of the 
constriction (see Fig. |l|) . The disorder strengths 7 corre- 
spond to mean free paths in the range 1.2 < kpl < 7000. 
The spatial distribution of impurities is the same for all 
curves. 




1 1.2 1.4 1.6 1.8 



L/L 

FIG. 4. Effects of disorder on the conductance of a wire 
of width 2k F R = 24 and initial length k F L = 60. The 
same configuration of impurities located both within the 
constriction and in the adjacent regions, as in Fig. |l[ is 
used for each curve. The disordered region is given by 
kpLdis = 36, and the impurity concentration is in all cases 
)<F n dis = 0.27. The impurity strengths are, from top to bot- 
tom, kpj = 0.29, 2.90, 7.26 and 21.8 ef, corresponding, re- 
spectively, to mean free paths k F £ = 7000,70,11 and 1.2. 
Curves are shifted by one unit of conductance. 

For a very large mean free path kpl = 7000 (solid line 
in Fig. the quantization is almost unaffected by the 
disorder. When I decreases, higher conductance plateaus 
are destroyed. Lower plateaus remain well defined at 
first, although they are shifted to lower values (see dot- 
ted line in Fig. [|). This shift is smallest for the lowest 
plateaus. For strong disorder, kpl ~ 1, the conductance 
quantization is no longer visible except for the very first 
step. 

ThiS|-|rbehavior is similar to some experimental 
results,0l23 where only a few conductance steps can be 
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observed, with a shift to conductance values lower than 
integer multiples of 2e 2 jh. This shift is experimentally 
observed to increase for higher conductance values, in 
agreement with our numerical results. 



found in experiments (see for example Refs. |q-|il|). Sub- 
tracting it from the total resistance, we indeed find that 
the peaks are shifted to quantized values, as is shown in 
the inset of Fig. [|. 



2. Conductance histograms 

As pointed out by several authors, EE conductance 
measurements are reproducible only in very controlled 
experiments and hence it is easier to study conductance 
quantization statistically in terms of histograms .B0O In 
Fig. H we show a conductance histogram representing 
380 different impurity configurations. The same geome- 
try, impurity concentration, and disorder strength as for 
the dotted line of Fig. || have been used, corresponding to 
a mean free path kpl = 70. Individual peaks (shown as 
dotted lines in Fig. ^| ) can be resolved by constructing 
histograms from partial conductance curves for elonga- 
tions L in the interval L n+ \ < L < L n , where L n is 
the midpoint of the n-th step of the clean system, i.e. 
G{L n ) = (2e 2 /h)(n-l/2). 




2 3 4 

Conductance [ 2e 2 /h ] 

FIG. 5. Conductance histogram made out of 380 individ- 
ual conductance curves. The width 27?, initial length Lo, 
length Ldis and impurity concentration m are chosen as in 
Fig. ^. The impurity strength is fcf.7 = 2.9 ef (mean free 
path kpl — 70). The inset shows the same histogram cor- 
rected by the calculated sheet resistance of R s ~ 585 fl. 

The peaks of the histogram are not located at integer 
multiples of 2e 2 /h, but are instead shifted to lower values, 
as is expected from the single curve of Fig. g. They can 
be moved back to quantized, values of the conductance, 
as for experimental results p by subtracting a classical 
resistance in series with that caused by the constriction. 
This additional contribution can be estimated using the 
Drude formula for the conductivity, a = , where n is 
the electronic density. For the mean free path chosen in 
the simulations we obtain a sheet resistance R s ~ 585 fl. 
This value of the resistance is consistent with the ones 
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FIG. 6. Mean values and widths of the peaks of the con- 
ductance histogram of Fig. |^. The bars show the numerical 
results, while the solid line gives the mean value as predicted 
from random matrix theory. The dotted line represents the 
conductance of the clean system and is shown for comparison. 

The resolution of individual peaks allows one to cal- 
culate the mean value G n and root-mean-square width 
AG n of the n-th peak. These quantities have been calcu- 
lated by Beenakker and Melsen using random-matrix the- 
ory (RMT) for a slightly different situation.E2l They have 
a set-up where (i) the disordered and constricted parts 
are spatially separated by scattering-free segments, (ii) 
the widening between the constricted and unconstricted 
parts occurs abruptly, and (Hi) the disorder is varied, 
while the constriction remains fixed. In our notation, 
their results for n open channels are given by 



riRMT _ AC 
u„ — — — 



2p 2 

AG« MT = 4- 



l+7r, 

IS' 1 



In 



1 + 7r, 



l + 6 7 „ 
(l + 7») f 



1/2 



(12) 
(13) 



where 7„ = (n + l)^f-i? s . These expressions are valid for 

^-Rs 1, i.e. for a good conductor. We can mimic the 
set-up of Ref. [2(] by sampling the conductance for elon- 
gations corresponding to the midpoints of the plateaus 
of the clean wire, where tunneling effects are negligible. 
The width AG l n is then exclusively due to the differ- 
ent impurity configurations. Fig. ^| shows our numerical 
results (bars), together with the RMT prediction (solid 
line), calculated from Eq. ( p^ ) using the same sheet re- 
sistance as in the simulations. The agreement between 
the mean values is excellent, and the widths, calculated 
from Eq. (|l3|), deviate less than one percent from the 
numerical results. This remarkable agreement may be 
due to the fact that the impurity concentration is low, so 
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that the constriction and the disordered parts are effec- 
tively separated by a scattering-free region, as assumed 
in Ref. 



20 



The hypothesis of abrupt widening between 
the constricted and unconstricted regions appears not to 
be crucial. 

TABLE I. Comparison of the different widths of the peaks 
defined in the text 



Peak # 


AG tun 




AG* MT 


AGlr + AG l „ 




1 


0.050 


0.102 


0.105 


0.152 


0.157 


2 


0.049 


0.149 


0.144 


0.198 


0.205 


3 


0.045 


0.177 


0.176 


0.222 


0.218 


4 


0.040 


0.203 


0.203 


0.243 


0.249 


•5 


0.035 


0.226 


0.225 


0.261 


0.257 


6 


0.028 


0.247 


0.244 


0.275 


0.276 



The histogram of Fig. || also contains a broadening 
AGj,*™ due to tunneling, which depends on the shape of 
the constriction and can be estimated from the results for 
the clean system. The total widths are found to be well 
reproduced by simple sums of impurity and tunneling 
contributions, AG* ot w AGj, + AG*™ , as shown in Table 

B . 

Finally, we want to emphasize that our result is clospt 
to experiment than previous numerical calculations, Ell 
in particular with respect to the peak heights, which 
decrease strongly with increasing conductance. We at- 
tribute the improved agreement to the fact that we con- 
sider an ensemble of different geometries, with fixed 
Fermi energy, as in the experiments, while in Ref. |2l] 
the Fermi energy was varied. Since .-these systems are 
nearly ballistic, the ergodic behaviorEl expected in the 
diffusive regime is violated: an ensemble of different ge- 
ometries with fixed cf is not statistically equivalent to 
an ensemble over cf with fixed geometry. 




L/L„ 



FIG. 7. Impurity effects on the tensile force for a wire of 
width 2&F-R = 12, initial length kpLo — 30 and a disorder 
characterized by kpLdis = 60 and kpi = 54 (upper and lower 
curves). The middle line represents the force of the corre- 
sponding clean system. 



3. Force 

In this last subsection, we study the effects of disor- 
der on the tensile force. In this case, it is no longer 
optimal to use a delta-function impurity potential. The 
force appears to be more sensitive to discretization than 
the conductance. In fact, the lattice constant a must be 
smaller than the range of the impurity potential. There- 
fore we use a Gaussian impurity potential with a finite 
extent d > a. 

The wire studied has a geometry given by kpR = 6 and 
kpLo — 30. The disordered length is kpLdis = 60 with 
an impurity concentration Xprii — 0.27, corresponding 
to a mean free path £ = 54. Results are shown in Fig. 
fj] for two different configurations of disorder. Although 
the disorder is rather weak, it has a strong effect on the 
tensile force, and can either weaken or strengthen the co- 
hesion. An intuitive understanding of these opposite be- 
haviors can be obtained by considering the change in the 
density of electrons due to disorder. An impurity located 
at the center of the constriction (cf. upper curve of Fig. 
0) would deplete the electron density in the constriction, 
and thus weaken the cohesion. Correspondingly, impu- 
rities located away from the center of the constriction 
tend to increase the density of electrons in the center, 
thus increasing the cohesion (cf. lowest curve of Fig. 
[?]). Further work is needed to understand these effects 
quantitatively. 



V. CONCLUSIONS 

In this paper, we have presented an exact numerical 
solution of the free-electron model of metallic nanowires,cl 
using recursive Green's function techniques. Our method 
allows, in principle, to calculate both the conductance 
and the cohesive properties for wires of arbitrary shape, 
but for computational reasons we have limited ourselves 
to two-dimensional wires containing a constriction with 
a simple geometry. 

The validity of the adiabatic and WKB approxima- 
tions used in Ref. [|was confirmed for long constrictions, 
while for short constrictions the cohesion was found to be 
enhanced. A persistence of the tensile force after closing 
of the last conducting channel has been observed. It is 
not yet clear if this is an important physical effect or just 
a curiosity due to the reduced dimensionality. 

We have commonly assumed an ad hoc shape r(z) of 
the constriction, which in general does not represent equi- 
librium. A more consistent way would be to calculate the 
equilibrium shape which minimizes the free energy for a 
given elongation. We have made a first step in this di- 
rection by minimizing the macroscopic part of the free 
energy. We have found that in this case the wire can be 
stretched more before it breaks. 

Extensive calculations have been made to investigate 
the effects of random impurities. Our model gives con- 
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ductance histograms very similar to those observed ex- 
perimentally. The peaks tend to decrease in height and 
broaden with increasing conductance. This confirms that 
disorder in the nanowire is the cause of the resistance that 
is usually subtracted from experimental results. The po- 
sitions and widths of the peaks are in excellent agreement 
with calculations based on random matrix theoryHa We 
have shown that the widths have two components, one 
due to impurities (perfectly described within random ma- 
trix theory), the other due to tunneling, which depends 
on the shape of the constriction. 

Our calculations have shown that the effect of disorder 
is particularly strong on cohesion, which is sensitive to 
the specific impurity distribution at the center of the con- 
striction. In fact, depending on the presence or absence 
of an impurity close to the center of the constriction, 
the cohesion is either decreased or increased. This find- 
ing, though qualitatively understandable, requires fur- 
ther study. 

Our model may appear to be oversimplified in several 
respects. We have used hard-wall boundary conditions, 
neglected both the electron-electron interactions and the 
ionic background, and assumed a specific geometry upon 
elongation. However, we do not think that these assump- 
tions are crucial for the essential effects studied in this pa- 
per. Other boundary conditions hardly affect the results. 
Coulomb interactions are found to be barely visible in the 
cohesionJ23L3-and are believed to preserve conductance 
quantizations As to the ionic background, one has to 
worry on the one hand, about its effects on the electronic 
structure, and on the other hand, about the kinetics of 
the deformation. The first effect ma*-,produce interest- 
ing fine structure in the conductance,!^ while the second 
can lead to-different paths for elongation and contraction 
processes,C-3 explaining the observed hysteresis.El A hys- 
teresis effect could be easily reproduced within our model 
by choosing different geometries for the two processes. 

Our approach, where the emphasis is on conduction 
electrons while the granularity of matter is neglected, 
is complementary to molecular dynamics simulations, 
where the atomic rearrangements are followed in detail, 
while the effects of conduction electrons are taken into 
account only p, an averaged way, e.g. by the embedded 
atom methodrJ The comparison between our results and 
experiments shows that, except for hysteresis effects, the 
overall behavior of both cohesive and electronic proper- 
ties of metallic nanoconstrictions are well captured by 
our very simple model. 
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